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We consider an extension of the Standard Model with an arbitrary number N of Higgs doublets 
(NHDM), and calculate their contribution to the oblique parameters S and T. We examine the pos¬ 
sible limitations on N from precision measurements of these parameters. In view of the complexity 
of the general case of NHDM, we analyze several benchmark scenarios for the Higgs mass spectrum, 
identifying the lightest CP-even Higgs with the Higgs-like particle recently observed at the LHC 
with the mass of ~ 125 GeV. The rest of the Higgses are put above the mass scale of ~ 600 GeV, 
below which the LHC experiments do not detect any Higgs-like signals except for the former famous 
one. We show that, in a scenario, with all the heavy Higgses degenerate at any scale, there are no 
limitations on the number N of the Higgs doublets. However, upper limits appear for certain not 
completely degenerate configurations of the heavy Higgses. 


I. INTRODUCTION 


The recent discovery of the ~H25 GeV scalar particle at the Large Hadron Collider (LHC) [TJ 2] perfectly fills the 
vacancy of the Higgs boson necessary for the completion of the Standard Model (SM) at the Fermi scale. Surprisingly, 
the SM with the Higgs boson in this mass range becomes formally self-consistent up to the Planck scale. In the 
absence of any signal of physics beyond the SM, this fact drastically strengthens the position of this model as the 
theoretical basis of particle physics. 

Although the new observed scalar state has so far all the properties expected of the SM Higgs boson, it is still 
possible that it could be a light scalar in a multi-Higgs extension of the SM, or a light supersymmetric Higgs boson, 
or a Higgs boson coming from a strongly interacting dynamics, where the theory becomes nonperturbative above the 
Fermi scale and the breaking is achieved through some condensate. Now the priority of the LHC experiments will be 
to measure precisely the couplings of the observed scalar to fermions and gauge bosons, and to establish its quantum 
numbers in order to identify it with one of these or some other options. On the other hand, searches for new particles 
beyond the SM are an essential task of the LHC experiments mm- 

In this paper, we consider a multi-Higgs extension of the SM, with an arbitrary number N of the Higgs electroweak 
doublets. Our goal is to study possible bounds on the number of Higgs doublets from the precision measurements of 
the oblique T and S parameters. 


We assume that the N Higgs 5(7(2) doublets are identical, with hypercharge equal to 1. Some features such as the 
relation between the mass and gauge eigenstates in the scalar sector and the relation of the Higgs vacuum expectation 
values with the symmetry breaking scale v ~ 246 GeV presented in the two Higgs doublet model are still fulfilled 
when the number of Higgs doublets is increased m- 


The paper is organize d as follows. In Sec. [TTJ we briefly describe the theoretical structure of the N Higgs doublet 
model (NHDM). In Sec. Ill we compute the one-loop contribution to the T and S parameters in the NHDM. The 
bounds on the number of Higgs doublets coming from T and S parameter constraints at 95%C.L. are computed in 
Sec. IV In Sec. [Vj we summarize our results. 


* Electronic address: antonio.carcamo@usm.cl 
tElectronic address: sergey. kovalenko@usm. cl 

* Electronic address: ivan.schmidt@usm.cl 





2 


II. THE MODEL 


We consider an extension of the SM with N copies of the complex SU(2)l weak doublet scalar Higgs fields with 
hypercharge Y = 1 (NHDM). The model scalar potential, invariant with respect to the SM gauge group, is 


N 




(!) 

ri 7 ^ 


i,j =1 


! N N 

+ 1 H E aijM ($ ir 2$.) (^r 2 ^) 1 

i,j,k,l=1 


i,j,k,l =1 


(i) 


where r 2 is a Pauli matrix in the SU (2) l space and 

&ij,kl — & ji,kl — (2) 

For simplicity, we assume all the parameters in the scalar potential to be real. Then the Hernriticity of the scalar 
potential 0 implies 

^ ij,kl ^ji,lkt ^ij,kl & ji,lk: Pij ft/'zi 

The minimum of the scalar potential is parametrized by N vacuum expectation values 


<*!> = 



We decompose the Higgs fields around this minimum as 
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( 5 ) 


where 


(ft) = (Vi) = (wj) = (■£*> = 0, / = 1,2,--- ,iV. (6) 

Then the covariant derivative acting on the Higgs doublets takes the form 

D^i = d,^, - l -gw;r a 

( + 775 b^ft - 0 W 2 («i + ft) + + g'YiB^) £ ; ] \ 
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2^2 V 2 - gW^ui + gWfa - (-<?W 3 + g 1 ^) {v t + Pl ) ) ’ 1 J 

where the r° are the ordinary SU(2)l Pauli matrices and V; = 1. 

The NHDM scalar-gauge boson interactions are given by 

E (^*0 = - E {2d M <ft + [ 5 ^ft - ( Vl + Pl ) + (<?W 3 + g'YtBj &] } 2 

1=1 0 1 = 1 

1 N 

+ g E i 2d ^i ~ {g w l + g' Y i B u) ui - [gWl {vi + Pi) + gWfa] } 2 

° 1=1 

1 N 

+ g { 2d vPi + [g w l£,i + gWfal - ( S W 3 - g'YiB^)] Pi} 2 
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1 N 

+ 8 + gWfa - (- 3 W 3 + g'YtBp) (v t + Pl )} 2 . ( 8 ) 


The connection between the interaction and mass scalar eigenstates is explained in what follows. The charged scalar 
fields of Eq. |5]) are linear combinations of the charged Goldstone bosons and the charged physical scalars. The 
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imaginary parts of the neutral component of the scalar doublets of Eq. ([5]) are linear combinations of the neutral 
Goldstone bosons and of the CP-odd neutral scalar fields. The real parts of the neutral component of the scalar 
doublets of Eq. ([5]) are linear combinations of the CP-odd neutral scalar fields. Within this framework we consider a 
scenario where the interaction and mass eigenstates are related in the way analogous to the two Higgs doublet model 
(2HDM) QIJ 


N 
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AT 
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( 11 ) 


Here u ss 246 GeV is the conventional electroweak symmetry breaking scale. The fields TV? (« = 1,2, ••• ,7V) and 
Aj (j = 1, 2, • • • , TV — 1) are the CP-even and CP-odd neutral Higgs bosons, respectively. Similarly to the W ± 
gauge bosons which are defined in terms of W 1 and W 2 , the charged Higgs and Goldstone bosons are related to the 
component fields in (10) as 


= ~ 3 


II) =F iH 2 


± 7T 1 =P ITT 2 


V2 ’ 


72 


( 12 ) 


Thus we assumed the following: 


1. The rotation matrix Q , which relates the neutral Goldstone boson n 0 and the CP odd neutral Higgses with the 
interaction eigenstate scalars (7 = 1, 2, • • • , TV) in Eq. is the same as the one that relates the components 


of the charged Goldstone bosons n 1 ’ 2 and Higgses Hj' 2 with the corresponding interaction eigenstates 


(7 = 1,2,-•• ,7V) in Eqs. (10), (12). 


2. The vacuum expectation values of TV Higgs fields vi (7 = 1,2, ■ 

breaking scale v ~ 246 GeV through the first relation in Eq. ©■ 


, TV) are related to the common symmetry 


Both assumptions are generalizations of the corresponding relations of the 2HDM m- In the case of NHDM, these 
relations are not true everywhere in the parametric space but only in a certain part of it. Adopting the above 
assumptions, we limit ourselves to a region in the parametric space of the NHDM, which is motivated (hinted) by the 
2HDM. 


III. ONE-LOOP CONTRIBUTION TO THE T AND S PARAMETERS. 


In this section we calculate one-loop contributions to the oblique parameters T and S defined as (T2UT7] : 


T = 


n 33 (q 2 ) - n n (q 2 




w 


s = 


2 sin 2d w dU 30 ( q 2 


q 2 = 0 


oiem{Mz) dq 2 


(13) 


q 2 — 0 


Here Hn (0), n 33 (0), and n 30 (g 2 ) are the vacuum polarization amplitudes with {WV W 3 }, { WW:) } and { W 2 , B^} 
external gauge bosons, respectively, where q is their momentum. Let us note that, in the aforementioned definitions 
of the oblique T and S parameters, it is assumed that the new physics is not light compared to M\y and Mz- 
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A. T parameter 


The interaction Lagrangian, relevant for the computation of one-loop contributions to the T parameter in Eq. (13), 
is 
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where P l: , is given by 
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By definition it satisfies the inequality 


0 < < 1 


(16) 


As seen from Eq. (141, the T parameter (13) at one-loop level receives contributions from the diagrams shown in 

Fig. [T] 

Their partial contributions, assuming the cutoff A to be much larger than the masses of the scalar particles, are 
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The subscripts in T 0 {, denote the internal lines of the diagrams in Fig.[T] The functions F (A 2 ,m 2 , m%) and G (A 2 , to 2 
are defined as 
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Figure 1: One-loop Feynman diagrams contributing to the T parameter. 


G (A 2 , to 2 ) 


lim F (A 2 , mi, m 2 ) = A 2 — 2 m 2 In 



+ m 2 


(24) 


Collecting all the contributions together, we find the one-loop contribution to the T parameter coming from the scalar 
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sector of the NHDM: 
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where we identified the lightest CP-even Higgs Pf = h with the LHC Higgs-like particle with the mass to?,, = 125 
GeV. 

The function /i (mf, m|) is given by: 


/i(m 2 ,m 2 ) = M^M 


ME, 


lim h (to 2 , to 2 ) = to 2 . 


(26) 


We can split the T parameter as T = Tsm + AT, where Tsm is the contribution from the SM, while AT contain all 
the contributions involving the heavy scalars: 
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B. 5 parameter 


The interaction Lagrangian relevant for the computation of the one-loop contribution to the S parameter in Eq. (13) 
is 
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As follows from this Lagrangian and the definition (13), the S parameter at one-loop level receives contributions from 
the diagrams shown in Fig. [2] 
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Figure 2: One-loop Feynman diagrams contributing to the S parameter. 


Their partial contributions, assuming the cutoff A to be much larger than the masses of the scalar particles, are 
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As before, the subscripts in S a b denote the internal lines of the diagrams in Fig. [2] Then, the 1-loop Higgs contribution 




























to the S parameter in the NHDM is 
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where we identified the lightest CP-even Higgs A? = h with the LHC Higgs-like particle with the mass nih = 125 
GeV. We defined a function 
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We can split the S' parameter as S = Ssm + AS, where Ssm is the contribution from the SM, while AS contain all 
the contributions involving the heavy scalars: 
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IV. T AND S BOUNDS ON NHDM 


The inclusion of the extra scalar particles modifies the SM predictions for the oblique parameters T and S, and 
therefore their values extracted from high precision measurements can be used to constrain the N Higgs extension 
of the SM. Our goal is to examine if these measurements are able to restrict the number N of Higgs electroweak 
doublets. The experimental results on T and S restrict the deviations AT and AS from the SM predictions to lie 
inside a region in the AS — AT plane. At the 95%C.L., this region is the elliptic contour shown in Fig. [3j taken 
from Ref. [TS]. The reference point AS = AT = 0 is conventionally taken to be the SM value of AS and AT at 
m,= 125.7GeV and m t = 173.18GeU. In view of the complexity of the general case of the N Higgs doublet model, 
we consider several benchmark scenarios described below. 


A. All the heavy Higgses are degenerate. 


This is the most simple case of the Higgs spectrum with the lightest CP-even Higgs A? = h identified with the LHC 
Higgs-like particle, with a mass mu = 125.7 GeV and all the other heavier Higgses degenerate having a common mass 
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Figure 3: The interior of the ellipse in AS — AT plane is the experimentally allowed region at 95%C.L. from Ref. [T8 . The 
reference point AS = AT = 0 is conventionally taken to the Standard Model value of AS and AT, at mh = 125.7GeV and 
mt = 173.18 GeV 


m#. Thus, 


m H ±=m A q=m H , j = 1, 2, • • • , N - 1, (39) 

3 3 

m H o = nih = 125.7 GeV, m H o = to#, i = 2, ■ ■ ■ ,N, 
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For this spectrum, Eqs. 


(28) and (38) for the AT and AS parameters 
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are drastically simplified and take the form 


fnl,m 2 H ) 


(40) 

(41) 


Thus, in this scenario neither of the two parameters S and T depends on N. Therefore, the spectrum in Eq. (39) 
does not constrain the number of Higgs doublets. 


B. Degeneracy inside the groups of the heavy CP-even, CP-odd and charged Higgses. 

Subscenario Bl: The CP-even and CP-odd neutral Higgses are degenerate. 

The next-to-simplest scenario that we consider has the following Higgs spectrum: 

m A o — m H , rn H ± = m H + A, j = 1, 2, • • • , TV - 1, (42) 

3 

m H o = mh = 125.7 GeV, m H o = to#, i = 2, • • ■ , N, 

m. H > m h - 


This spectrum, using Eqs. (281 and (38), leads to the expressions 
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Figure 4: Upper bound N < lV max on the number N of Higgs doublets, obtained from T and S at 95%C.L., using the 
experimental constraints indicated in Fig. [3] for the Higgs spectrum in Eq. (|42[). 
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where we used the first relation in Eq. @. 

Now, using Eqs. (43), we find the maximal values N n 
in Figi We scan the parameter space within 


of the Higgs doublets N compatible with the precision data 


0 < Pn < 1, 600GeV <m H < ITeV. 


(45) 


In Fig. [4j we show the resulting values N max in a function of the splitting parameter A. It is noteworthy that 
the maximum number of Higgs doublets decreases when the mass splitting A between the heavy physical scalars is 
increased. This behavior follows from the fact that increasing the number of Higgs doublets yields an increase of the 
T and S oblique parameters. 


In the limit A —> 0, we find N rnax -A oo, corresponding to no limits on N, which is consistence with the scenario (391. 
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Subscenario B2: The CP-even neutral and charged Higgses are degenerate: 

m H ±=m H , m A o=m H + A, j = 1,2, • • • , N - 1, (46) 

j 

m H o = mh = 125.7 GeV, rn H o = ran, i = 2, • • ■ ,N, 

m H > m h - 


N 

' v max 



Figure 5: The same as in Fig. [4] but for the spectrum in Eq. (|46|). 


From Eqs. (28) and (|38|) , we find for this spectrum 
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Scanning the parameter space in the region (|45j), we find the maximal values on the number N of Higgs doublets 
compatible with the data in Fig. [3j The results are shown in Fig. [5| As seen, this spectrum is significantly less 
restrictive for N than that in Eq. (42). This is mainly because of the fact that only the S parameter depends on N 
in the present case, while in the case of the spectrum (42) both T and S are A^-dependent. 


Subscenario B3: The CP-odd neutral and charged Higgses are degenerate: 
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In this case we find from Eqs. 
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Figure 6: The same as in Fig. [4] 
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values on the number N of the Higgs doublets compatible with the data of Fig. [3] The results are shown in Fig. [6] 
Again, as in the case (461, the parameter T is independent of N. As a consequence, the limits on N for the spectrum 


(48) are significantly weaker than for (42). 


Subscenario B4: Split groups with the degenerate interior: 
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For all these cases we find from Eqs. (28) and (38) 
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As seen from Eq. ( |53| , now both parameters T and S depend linearly on TV, the number of Higgs doublets. Scanning 
the parameter space in the region (45), we find the results for several sample values of the parameter a = 0.5,1,2 
shown in Figs. m Note the general tendency: the large splitting, corresponding to the larger values of A and a, 
leads to more stringent constraints on TV. 



Figure 7: The same as in Fig. |4j but for the spectrum 
(dashed line), a = 2 (solid line). 


in Eq. (501. 


The left panel is for a 


1; the right panel is for a = 0.5 


C. No degeneracy with a particular structure of the spectrum 


Finally, let us consider a benchmark scenario in which all the Higgses are nondegenerate. Since the general case can 
hardly be analyzed, we consider a particular structure of the Higgs spectrum: 

m H° = mh = 125.7GeE, m H o = + 2 (i — 2) A, i = 2, • • ■ , TV, (54) 

m A o = m// + (2j — 1) A, m H ± = m H + 2 (j - 1) A T S, j = 1, 2, • • • , TV - 1, 

J j 

777TT > 777/j. 

This is a spectrum equidistant in each of the three groups of Higgses {A^}, {T?°}, {Hf} with a step 2A. The bands 
of these groups overlap with each other. The spectrum is characterized by two parameters A and <5. 
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^max ^max 



A(GeV) 


Figure 8 : The same as in Fig. [?] but for the spectrum in Eq. ( |51| ) . 


For simplicity, we assume 


P ■ — S- 

r i,3 ~ u *3i 


i,j = * = 1,2,--- ,N 


(55) 


Then from Eqs. (J28J) and (|38j) we get 
AT ~ 


N-l 


16ir 2 v 2 a E M{M z ) 


T, \rn 2 H ± - h (m 2 A o,rn 2 H ±^J + h (m 2 H o +i ,rn A o) - h ( m ff° +1 i m jy±) 


N -1 


as 1 ~ —h y K 

127T 


i=l 
, TTlj^o , 


)■ 


(56) 

(57) 


Scanning the parameter space in the region (451, we find the results for several sample values of the parameter S = A/n 
shown in Fig. [To] The curves from the bottom to the top correspond to n = 20,10, 2, respectively. With the larger 
value n = 50, we find N max ~ 570 for A = 20 GeV. 


V. SUMMARY AND CONCLUSIONS 


We have considered an N Higgs SU( 2) doublet model (NHDM) with arbitrary number N. In this model, we calculated 
the one-loop contributions AS and AT of the Higgs doublets to the electroweak oblique parameters S and T. The 
calculated contribution depends on the number N of Higgs doublets, and, therefore, our results can be used to 
constrain N from data on the precision measurements of the parameters T and S. Within the generic case of the 
NHDM, due to the large number of free parameters, this program can hardly be realized. For this reason, we have 
analyzed several benchmark scenarios with particular mass spectra [Eqs. 


(42), (46), (48), 


([50j)-(|52]), and Q] of 
mown case of 


the physical scalars of the NHDM, including some other simplifying assumptions, inspired by the well-, 
the 2HDM, about the physical Higgs mixing and the vacuum structure of the model. These scenarios correspond to 
certain domains of the NHDM parameter space. We have shown that, except for a very particular “fine-tuned” case 
with all the physical heavy Higgses degenerate [Eq. (39)], these scenarios imply constraints on the number of Higgs 
doublets N 1 in order to be compatible with the existing data on the precision measurements of T and S. 
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^max 



Figure 9: The same as in 
a = 0.5,1, 2, respectively. 


Fig. a but for the spectrum in Eq. 


The curves from the bottom to the top correspond to 


We presented our results on TV < N max in Figs. m as functions of the mass splitting parameter A. The general 
feature of our results is that the maximal number N max of Higgs doublets is a monotonically increasing function for 
small values of the splitting A < 20 GeV and a monotonically decreasing one for larger values. Thus, the data on 
T and S are able to accommodate an arbitrary large number TV with decreasing splitting between the masses of the 
physical scalars, and, vice versa, TV becomes stringently constrained in those parts of the NHDM parameter space 
with large mass splitting in the scalar sector. The same tendency is demonstrated by the plots with respect to their 
dependence on the two other additional parameters a and <5, characterizing the Higgs mass spectrum: a smaller mass 
splitting corresponds to a larger number TV compatible with the analyzed data and visa versa. This is the main 
message of the present study. Also worth mentioning is that the maximal number N max of Higgs doublets is exactly 
the same when the charged Higgses are degenerate between either the CP-even or the CP-odd neutral Higgses, as 
shown in Figs. |5|6[ and larger than the obtained in the scenario of CP-odd - CP-even neutral Higgses degeneracy. 
Consequently, the tightest T and S oblique parameter constraints arise in the scenario where the charged Higgses 
are split in mass between either of the neutral CP-even or CP-odd Higgses, as shown in Fig. [4] This indicates that 
making the charged Higgses degenerate between either the neutral CP-even or the neutral CP odd Higgses helps avoid 
constraints from electroweak precision observables, a feature already present in the 2HDM. 


Our analysis cannot exclude a deviation from this tendency in certain parts of the NHDM parameter space; however, 
in our opinion, this should be related with certain “fine-tuning” of the parameters, as in the case of the spectrum 
(39). Naturally, the concrete limit on TV depends on a particular scenario within the generic NHDM framework. We 
hope our results will help examine such scenarios regarding their consistence with the present and future data on the 
precision measurements of the electroweak oblique parameters T and S. 


As a final remark, let us note that the maximum number of Higgs multiplets is constrained from the requirement 
of perturbative unitarity in the scattering amplitudes of two transversely polarized W bosons into a scalar pair, as 
shown in detail in Ref. m From the results given in Ref. [19] , it follows that the maximum number of Higgs 
doublets consistent with the aforementioned perturbative unitarity requirement is 2307, which is larger than our 
obtained upper bounds on the number of Higgs doublets, from oblique T and S parameter constraints, for the several 
benchmark scenarios considered in this paper. Consequently, our upper bounds on the number of Higgs doublets are 
consistent with the perturbative unitarity in the scattering amplitudes of two transversely polarized W bosons into a 
scalar pair. 
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Figure 10: The same as in 
n = 20,10, 2, respectively. 


Fig. a but for the spectrum Eq. (|54[). 


The curves from the bottom to the top correspond to 
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